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Abstract — In this paper, we propose a novel adaptive modulation 
and coding (AMC) algorithm dedicated to reduce the feedback 
frequency of the channel state information (CSI). There have been 
already plenty of works on AMC so as to exploit the bandwidth 
more efficiently with the CSI feedback to the transmitter. However, 
in some occasions, frequent CSI feedback is not favorable in 
these systems. This work considers finite-state Markov chain 
(FSMC) based channel prediction to alleviate the feedback while 
maximizing the overall throughput. We derive the close-form of the 
frame error rate (FER) based on channel prediction using limited 
CSI feedback. In addition, instead of switching settings according 
to the CSI, we also provide means to combine both CSI and 
FER as the switching parameter. Numerical results illustrate that 
the average throughput of the proposed algorithm has significant 
performance improvement over fixed modulation and coding while 
the CSI feedback being largely reduced. 

I. Introduction 

In recent years, the demand of wireless communication 
channel has grown tremendously. It may be not sufficient to 
support the large amount of data over a bandwidth-limited time- 
varying fading channel using fixed modulation and coding. The 
fixed modulation and coding schemes with power control are 
designed to maintain certain quality of service (QoS) when 
the channel condition is poor. Thus, fixed modulation and 
coding scheme does not fully exploit the time-varying nature 
of the channel. Adaptive modulation and coding (AMC), on the 
other hand, takes advantage of the time-varying characteristic 
to achieve high spectral efficiency while maintaining acceptable 
bit error rate (BER). The notion of adaptive transmission has 
been proposed striving to satisfy growing consumer demands 
(TJ. Moreover, some adaptive transmission techniques have been 
widely included in wireless standards such as IEEE 802.16 Q, 
IEEE 802.22 Q, and LTE g). 

Several pioneer works on AMC have been proposed in the 
literature. In J5j and (6), various rate and power adaptation 
is considered. Goldsmith et. al. 121 proposes a variable-power 
variable-rate modulation scheme. The transmitter switches set- 
tings according to the channel state information (CSI) feedback 
by the receiver. The transmitter compares the received channel 
gain to a set of prescribed thresholds. The thresholds partition 
the channel gain into a set of regions corresponding to dif- 
ferent adaptation settings, and the transmitter makes decision 
accordingly. Once a setting has been determined, the setting 
is unchanged for a transmission frame. The settings in J7) 
and JSj are restricted to variable-rate while the power is kept 
constant. In [9|, a frame consists of a variable number of packets 
depending on the current channel state. The settings of ||7ll-ll9l 



are switched at the start of transmission of each frame. All 
of these works require current channel state information at the 
receiver, thus a reliable feedback path is required for the CSI 
feedback. Moreover, the receiver is obliged to feedback the CSI 
every frame. 

However, in some circumstances, it is prohibitive to maintain 
a high feedback rate. For the systems where power is a concern, 
such as sensor networks [ 1 1 , frequently CSI feedback is not 
available. Moreover, the research on adaptive transmission tech- 
nique combined with multiple-input multiple-output (MIMO) 
or orthogonal frequency division multiplexing (OFDM) or both 
of them has already received more and more attention. The 
adaptive OFDM system [ 1 1 J — [ 1 3 1 involves optimizing the mod- 
ulation level and transmission power over the entire frequency 
band, so the CSI is required for all subcarriers resulting in high 
feedback load. In MIMO systems lfT4l . the multiple antennas 
lead to the increasing complexity of CSI. As a result, the CSI 
has higher dimension and is more complex in these systems 
than that in the single carrier case, and there have been efforts 
on the reduction of feedback load in these system. 

Therefore, the further reduction of feedback overhead is 
desired. One approach is to reduce the frequency of CSI 
feedback from every frame to every M frames. The adaptation 
period of the system is also reduced to every M frames. 
However, in the previous works, adaptive transmission makes 
decision based on current channel state [5|-[8|, or the most 
probable channel state estimated by ACK/NAK signal ifTSl . If 
the adaptation period is reduced to every M frames, the number 
of the CSI received at the transmitter side decreases. Without 
the knowledge of channel variation between intervals of CSI 
feedback, the adaptive method is chosen only to maximize 
the throughput of the first frame of the M frames instead of 
the throughput of all subsequent M frames. Consequently, the 
performance of such system is poorer than the system with 
shorter adaptation period. 

In order to maximize the overall throughput, it is required 
to take the channel state of subsequent M — 1 frames into 
consideration. Thus, some predictions on the channel state of 
subsequent M — 1 frames are to be made. Currently there are 
several channel prediction method such as linear prediction 1 16 1, 
ifTTl . channel variation prediction based on finite state Markov 
chain (FSMC) flj], ED, etc. However, in the work lfl3l. Ifl8l. 
only the channel state prediction of the next frame is considered. 
The settings are selected to maximize certain criteria of only the 
transmission of the next frame. Their results could not extend 
the adaptation period to M frames. 



The previous works on channel prediction does not provide 
the analysis of BER or FER over long variation of channel. In 
this paper, a Markov chain-based channel prediction method is 
proposed. We are able to predict the channel variation of the 
following M — 1 frames through the first frame and derive the 
close-form expression of the frame error probability of subse- 
quent M — 1 frames. Hence the channel variation of subsequent 
M — 1 frames are taken into consideration, which enables the 
optimization the overall throughput of the M frames. On the 
basis of channel prediction based on FSMC, we propose a novel 
adaptation algorithm: Prediction-based Adaptation (PRADA) 
Algorithm for modulation and coding. The proposed algorithm 
is able to extend the adaptation period from every frame to 
every M frames while still maintaining good performance. 
The numerical result shows that with channel prediction, the 
throughput of system using PRADA increases significantly over 
the system that only maximizes the throughput of the first frame. 

The systems proposed by the previous works 0-||9) switch 
settings based on the received channel gain. There are also 
works proposing to utilize solely ACK and NAK feedback 
03]. Another approach does not assume the channel condition 
but simply use ACK and NAK as setting switching criterion 
020 ■ The transmitter first selects the largest constellation which 
has the highest data rate and the highest error probability, 
and reduces the constellation size once receives retransmission 
request. The another approach uses ACK and NAK as an ob- 
servation to estimate the underlying channel state lfT31 . Il20l . In 
ITT31 the transmitter forms a belief state space using underlying 
FSMC model. The transmitter determines the most possible state 
according to the received ACKs or NAKs, and switches settings 
accordingly. On the other hand, in the algorithm proposed by 
Vucetic ll20l . the transmitter compares the block errors to the 
specified thresholds to estimate the current channel state. These 
thresholds are determined by minimizing the errorneous channel 
estimation. Both algorithms adopted the FEC codes of different 
rates correspond to the estimated channel state. 

In this paper, we develop another adaptation algorithm 
PRADA-B which combines both feedback-reduced CSI and 
ACK7NAK signal, while the precedent works use either CSI or 
ACK/NAK only. The transmitter compares the number of frame 
errors in M frames to the thresholds which are computed by 
the CSI feedback every KM frames. The thresholds are deter- 
mined with the view to maximize the overall throughput over 
consecutive KM frames, instead of optimization on the current 
frame. Switching settings according to FER can further reduce 
the CSI feedback while maintaining acceptable throughput. The 
receiver feeds back the CSI to the transmitter at the beginning 
of every KM frames in order to track the channel behavior. 

The rest of this paper is organized as follows. The system and 
channel model is given in Section II. The adaptation algorithm 
of PRADA is given in Section III. We develop the mathematical 
basis of the algorithm, the frame error rate derivation and the 
threshold of setting switching in Section IV. The numerical 
results of the proposed algorithm are demonstrated in Section 
V. Finally, we draw concluding remarks in Section VI. 
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II. System Model 

A. System Overview 

The system model considered in this work is given in Figure 
Q] The information bits are encoded by the rate-compatible 
punctured convolutional (RCPC) code, modulated, then trans- 
mitted. The transmitter switches the modulation and the coding 
every setting adaptation period. The setting adaptation period is 
M frame periods. 

In this model, R settings are chosen from the combina- 
tions of modulations and RCPC codes. Denote R settings as 
{s%, S2, ■ ■ ■ , sr}, listed in the descending order of the spectral 
efficiency, which has the same meaning as listing the settings 
in the descending order of FER in this system. The FER is 
computed based on the same number of modulation symbols. 
Each setting corresponds to a data rate, which combines the 
effect of the modulation and the coding. Denote the data rate of 
each setting as {fci, kz, . . . , k^}. Note that k\ > k^ > ■ ■ ■ > &r 
because the setting with higher data rate also has higher FER. 

The receiver sends an ACK or a NAK to the transmitter when 
it receives a frame. In addition, the receiver sends the CSI to 
the transmitter at certain frequency. The transmitter switches the 
setting according to the CSI and the FER. The algorithms are 
described in Section ITlTl 

B. Channel Model 

FSMC is widely used to model Rayleigh fading channel 
which is a good approximation of channel with memory. And in 
the application of packet/frame transmission, first-order Markov 
chain's accuracy is supported by [21 1. There are some studies 
on how to partition SNR into finite states of FSMC model. 
Wang and Moayeri [22 1 propose a method partitioning all 
FSMC states that have the same stationary probability. Zang 
and Kassam [23 1 propose a method partitioning all FSMC states 
that have the same average time duration. Park and Hwang 



propose a method partitioning FSMC states such that the 
non-adjacent transition probabilities are small enough. In this 
work, we adopt the channel model proposed in |f23l which is a 
better approximation and the average time duration is made large 
enough to meet our assumption that in each frame transmission 
the state will not change. 

The channel SNR 7 is quantized into several non-overlapping 
states with boundary values T — [Til^ . . . Tjv+i] with Ti = 
and Tat +1 = 00. The channel is said to be in state Wi if 
7 G [r,Tj_|_i). Thus we can obtain stationary probability 7T, of 
Wi using the probability density function of Rayleigh fading 
channel. 

f Tt+1 1 7 

7Ti = / — exp( )d7 (1) 

JFt 7o 7o 

where 70 is the average SNR. 

In this system, we assume constant time in each frame tran- 
sition, denoted by T p . We can obtain the transition probabilities 



M frames 



(2) 



where AT(I\-) is the level crossing rate at SNR level I\ either 
in the positive direction or in the negative direction, given by 



N(T) = J—f m ex V (--) 

V 70 70 



(3) 



where f m denotes the maximum Doppler frequency due to the 
user motion. 

C. FER Model 

Given the channel model, we can analyze the performance of 
the system. Let P Sr {pi) denote the FER of the channel SNR 7 
when using setting s r . The FER of the channel state Wi when 
using setting s r , which is denoted by Ef\ can be derived as 



(4) 



where ptil) is the probability that channel SNR is 7. 

In this work, the FER is obtained by computer simulation 
using RCPC codes. The modulation and the puncturing rate 
are corresponding to the setting s r . However, the expression 
can also be computed using analytical expression or simulation- 
fitting expression of P Sr {l)- 

III. Algorithm Descriptions 

In this section, we describe the proposed algorithms, that 
is, PRADA-A and PRADA-B. PRADA-B considers the CSI 
feedback and the FER in previous transmissions. In order to 
choose the optimal setting, it requires a threshold table for FER 
comparison. Although the computation can be done offline, the 
complexity may be not allowed in some system. PRADA-A is 
also proposed because it reduces the computational complexity 
of PRADA-B. PRADA-A considers only the CSI feedback and 
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Fig. 2. Time domain diagram of Prada-A algorithm 
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Fig. 3. Time domain diagram of Prada-B algorithm 



its computation is only the first part of PRADA-B 's computa- 
tion. 

By simulation, we can get the FER of each setting in each 
state. The expected throughput of successive M frames given 
the first channel state and the setting can be derived if the FER 
and the channel transition probability are known. It is used in 
both PRADA-A and PRADA-B. The derivation of the expected 
throughput is given in Section HV-BI 

To deal with the doppler frequency variation, the algorithms 
use different FSMC channels. That is, we can generate several 
FSMC channels for possible doppler frequency and use them in 
algorithm computation. 

A. PRADA-A Algorithm 

Figure [2] illustrates the time domain diagram of PRADA-A. 
We call M frame transmission duration a setting adaptation 
period. At the beginning of a setting adaptation period, the 
transmitter chooses the setting maximizing the throughput of 
the next M frames according to the current channel state. Here 
we summarize the steps as follows: 

• Step 1 : For every setting, compute the expected throughput 
of the next M frame transmissions for every possible 
beginning channel state. This can be done offline. 

• Step 2: At the beginning of a setting adaptation period, 
the transmitter observes the channel state and chooses the 
corresponding setting. 

• Step 3: Repeat 2. 

B. PRADA-B Algorithm 

Since the FER in the previous setting adaptation period 
implies the wireless channel condition implicitly, we propose 
a FER-based algorithms. 

The time domain diagram of this algorithm is shown in 
Figure [3] A threshold adaptation period is K setting adaptation 
periods. At the end of a setting adaptation period, the transmitter 



compares the FER to some predetermined threshold to decide 
which setting to used in next setting adaptation period. At 
the beginning of a threshold adaptation period, the transmitter 
decide the threshold used in the next K setting adaptation 
periods. The detail of the threshold determination is given in 
Section |IV] 

The adaptation consists of two parts, one is the adaptive 
modulation and coding switching, and the other is the threshold 
adaptation. The former is for the short-term channel variation 
while the latter is for the long-term channel variation. 

1 ) Adaptive Modulation and Coding Switching: The trans- 
mitter computes the FER in the previous setting adaptation 
period first, and then compare the FER to the threshold. A 
threshold set consists of 2R values {T r r -i,T r . r+ i}^ =1 , where 
Ti, = and T R<R+ i = 1. If the FER is lower than T rjr _i 
when using setting s r , the transmitter switches to setting s r _i; 
if the FER is higher than T r r +i when using setting s r , the 
transmitter switches to setting s r +i; otherwise, the transmitter 
uses the same setting as before. At the beginning, the system can 
start with any setting because it would switches to the optimal 
setting later. 

The steps are summarized as follows: 

• Step 1: Initiate the transmitter using setting s%. 

• Step 2: In each setting adaptation period, record number 
of erroneous frames I, and the FER = l/M. 

• Step 3:If the FER is lower than T r r _ 1 , switch to setting 
s r _i; else if the FER is higher than T r r+ i, switch to 
setting s r+ i\ else the setting remains the same. 

2) Threshold Adaptation: Due to the large-scale fading of 
wireless channel, the transmitter changes the thresholds every 
K setting adaptation periods (i.e., KM frame transmissions). 
The thresholds are chosen to maximize the expected throughput 
of the next KM frame transmissions. The transmitter changes 
the thresholds at the beginning of the KM frame transmissions 
based on the channel state Wi and the current setting s r . The 
detail of computing the thresholds is given in Section [IV] 

Here we summarize the steps as follows: 

• Step 1: For each possible combination {s r Wi}, find the 
corresponding threshold set. This can be done offline. 

• Step 2: At the beginning, use the threshold set correspond- 
ing to {siWi} for the first KM frame transmissions, where 
Wi is the channel state at the beginning. 

• Step 3: After K setting adaptation periods, the system 
changes the threshold set according to the first state {s r Wi} 
of the next K setting adaptation period. 

> Step 4: Repeat 3. 

IV. Algorithm Analysis and Derivation 

In this section, we describe the analysis of the PRADA 
algorithm. Since PRADA-A is a simplified version of PRADA- 
B, it only related to Section IIV-BI 

A. Induced Markov Chain 

Let P w denote the channel-state transition matrix, where P/" 
is the probability that channel state changes from i to j in a 
frame transmission. Let F™ k (s r ,l) be the probability that in a 
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switching 



setting adaptation period, channel state ends at Wk and I frames 
are erroneous when using setting s r given that channel state 
starts from Wi. The derivation of F™ k (s r , I) is given in Appendix 

Let {A t } denote the induced Markov chain which models 
both setting switching and the varying wireless channel states. 
Each state represents the first state of a setting adaptation 
period. Denote the state space of {A t } by 5a = {s r Wi\r € 

{1,2 i2},< G {1,2, ... , N}} = {A (m) } where A (m) = 

s (m.) w (m) denotes the m-th element of the indexed state space. 
Here we derive the transition probability matrix of A t . 

If the first state of this setting adaptation period is s r Wi and 
the first state of the next setting adaptation period is s r _iiUj-, it 
means: 

1) the number of erroneous frames I in this setting adaptation 



period satisfies that l/M < T r r _i 
2) the channel state changes from some state to Wj at the 
beginning of the next setting adaptation period. 

Hence, the transition probability can be derived as 

k=l 1=0 

We can derive P SrWi ,s r w and Ps rWi ,s r+1 w by the same method. 
Thus, the transition probabilities can be express as 

N |MT,,,._iJ 

Ps r Wi,s r -iWj — ^ ' ^ ] Fi,k( s r,l)Pk,j 
k=l 1=0 
N |M2V,r+lJ 

P SrWi , SrWj =J2 E F Usr,l)P^ (5) 

fc=i ;=LMT rjr _ij+i 

N M 

Ps r Wi,s r+1 w, = y ' ^ ] Pi.k{ s rT^)Pk,j 

k=l l=[MT rtr+1 \ + l 

As a result, we get the transition probability matrix of the 
induced Markov chain P. The transition diagram are shown 
in Figure [4] 

B. Expected Throughput 

Let £, Sr wi denote the expected throughput of M frame trans- 
missions when the first state is s r Wi. Recall that R settings 
{si, S2, ■ ■ ■ , sr} are in the descending order of the FER. Each 
setting s r corresponds to a data rate k r which combines the rate 
of modulation and coding. In general, {k\, k2, ■ ■ ■ , k^} are also 
in the descending order. For each frame, we define its expected 
throughput as k r (l — Ef^. Hence, the expected throughput for 
M frames is 



N M M-l 

£,s rWi — k r — F™ k (s r ,i). 

k=l 1=0 



(6) 



Note that F™ k (s ri l) contains the probability that the channel 
state ends at w^. and the probability that I frames are in error. 



C. Threshold 

The throughput of KM frames is defined as 

K 



i K 



(7) 
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where ^ is the average throughput of i-th period. The expected 
throughput of KM frames given Ai is 
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0.9736 
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4 


0.0177 


0.9737 
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5 


0.0187 


0.9737 


0.0076 
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0.0198 


0.9736 


0.0066 
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0.0159 


0.9841 


0.0000 



TABLE II 

Transition probabilities of the FSMC. 
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2.0499 
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5.6514 
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7.0454 
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TABLE III 
SNR THRESHOLD VALUES OF THE FSMC 
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\(j-)\k m = A(j)] is the n-step 



where 

A(i),A(j) 

transition probability and 4A (m) is given in Equation (|6). Recall 
that P^™- 1 = P" only depends on the threshold set, as that 
described in Section HV-AI Thus, we can choose the threshold 
set such that P[£t|Ai] is maximized for every possible Ai. 

An algorithm used to find the threshold set is described as 
following. For each possible Ai 

1) Generate a threshold set {T,. ir _i, T Ttr+ i}^ =1 as the cur- 
rent threshold set randomly. 

2) Change one value in the current threshold set by 1 or —1 
to get a new threshold set. Do this for each value except 
for Ti q and Tr and get 2R — 2 new threshold sets. 

3) Compute the expect throughput for 2R — 2 new threshold 
sets. Choose the threshold set that maximizes the expect 
throughput to be the current threshold set. 

4) Repeat step 2-3 until all 2R — 2 new threshold sets are 
not better than the current threshold set. Then output the 
current threshold set. 

5) Do step 1-4 several times, choose the best threshold set 
among the output threshold sets. 
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Modulation 


Code rate 


Frame size 


Data bits/Frame (k r ) 


si 


16-QAM 


2/3 


8192 


5461 


S2 


16-QAM 


1/2 


8192 


4096 


S3 


QPSK 


2/3 


4096 


2731 


S4 


QPSK 


1/2 


4096 


2048 


S5 


BPSK 


1/2 


2048 


1024 



TABLE IV 

THE SETTINGS USED IN THE SIMULATION. THE FRAME SIZE CORRESPONDS 
TO 2048 SYMBOLS. 



t=l m=l 



State si S2 S3 S4 S5 



1 


1.0000 


1.0000 


1.0000 


1.0000 


1 .0000 


2 


1.0000 


1.0000 


0.7741 


0.5406 


0.3662 


3 


1.0000 


0.9387 


0.0900 


0.0229 


0.0128 
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0.9193 


0.3752 


0.0020 


0.0000 


0.0000 


5 


0.3807 


0.0411 


0.0000 


0.0000 


0.0000 


6 


0.0609 


0.0026 


0.0000 


0.0000 


0.0000 


7 


0.0045 


0.0000 


0.0000 


0.0000 


0.0000 



The zero values are actually very small numbers thus 
neglected. 



TABLE V 

FER OF EACH SETTING IN DIFFERENT CHANNEL STATES 



Scenario 


Average throughput 
(Data bits/Frame) 


PRADA-A 


1829.7 


PRADA-B 


1873.9 


16-QAM (2/3) 


527.6 


16-QAM (1/2) 


920.8 


QPSK (2/3) 


1724.8 


QPSK (1/2) 


1544.4 


BPSK (1/2) 


854.2 



TABLE VI 

The average throughputs of different scenarios. 



V. Numerical Result 

A. Performance of Different Average SNR 

In this channel model, the doppler frequency is set to 4 Hz 
and the average SNR is set to 2 dB. The FSMC used to simulate 
the wireless channel is shown in Table [TTJ The SNR boundary 
values in the FSMC are computed according to Section III-BI 
listed in Table HTT1 where state Wi is between Tj and Ti+\. When 
channel is in w\, the transmitter does not transmit signals. 

In the simulation, we consider five settings. Table [IV] shows 
the details of the five settings. Note that the frame durations of 
all settings are equal. The simulated FER values are listed in 
Table El 

To evaluate the performance of the PRADA algorithm, we 
simulate the adaptation case (PRADA-A, M = 120; PRADA- 

B, K = 4, M = 30) and the no adaptation cases (si, S2, 
s 3 , and s 5 ). We set the PRADA-A and the PRADA-B to 
have the equal frequency to obtain channel state information. 
That is, the feedback rates of both PRADA-A and PRADA- 
B are the same (per 120 frames). For fair comparisons, all 
simulation scenarios experience the same channel variations. In 
the simulation, 5,000,000 frames are transmitted, the average 
throughputs over time are shown in Table [VI] 

The result shows that the performance of the PRADA is better 
than that of the fixed modulation. The throughput over time is 
shown in Figure [5] We observe that while the performance of 
the fixed modulation saturates, PRADA algorithm can switch to 
the high-rate setting. When the channel is in bad condition, the 
PRADA algorithm can switch to the low-rate setting to maintain 
an acceptable performance. 

Figure [6] shows the throughput cumulative distribution func- 
tion (CDF) comparison of all cases. The CDF curves of the 
PRADA is better than other curves. At the low throughput, the 
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Fig. 5. Throughput over time. Each point is the average throughput of 30 
frames. 
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Fig. 6. CDF comparison of all cases 

curve of the PRADA is close to the curve of the QPSK (1/2); 
at the median throughput, the curve of the PRADA is close to 
the curve of QPSK (2/3); at the high performance, the curve of 
the PRADA is close to the curve of 16QAM (2/3). This implies 
the adaptation of the PRADA algorithm. 

Throughputs of different average SNR are plotted in Figure 
[7] To get the gain of PRADA algorithm, we compute the ratio 
of the throughput of PRADA to the highest throughput among 
the five fixed settings. Figure [8] shows the throughput gain of 
PRADA-A and PRADA-B. 

B. Performance of Different Doppler Frequencies 

In this simulation, the average channel SNR is set to 2dB, and 
the doppler frequency is from 2 Hz to 20 Hz. This simulation 
compares the following three algorithms. 

1) PRADA-A (M = 120) 

2) PRADA-B (M = 30, K = 4) 
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Fig. 7. Throughputs of different average SNR 
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Throughput (data bits/frame) 


PRADA-A 


1768 


PRADA-B 


1792 


Compared Algorithm 


1731 



TABLE VIE 

The throughput with doppler frequency variation 
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Fig. 8. Throughput gain of PRADA algorithm 



State index 


pw 

r i,i-l 


pw 

M+l 


1 





0.00268/ 


2 


0.00386/ 


0.00271/ 


3 


0.00416/ 


0.00242/ 


4 


0.00442/ 


0.00216/ 


5 


0.00468/ 


0.00190/ 


6 


0.00493/ 


0.00164/ 


7 


0.00396/ 






/ is the doppler frequency in Hz. 



1-(P« 



-1 



hi) 



TABLE VII 
Transition probabilities of the FSMC. 



3) A compared AMC algorithm that chooses the setting such 
that the throughput of the first frame is maximized. The 
setting is switched every 120 frames. 
Note that three algorithms above have the same rate of the CSI 
feedback. The settings used in the simulation is listed in Table 
HVl The FSMC channel boundaries are listed in Table [TTT] and 
the transition probabilities are list in Tabl dVIll 

Figure [8] shows the results of different doppler frequencies. 
PRADA-B is always better than other two algorithm. When 
the doppler frequency becomes high, the performances of three 
algorithms decrease. However, PRADA-A and PRADA-B at 
least maintain a performance that the best fixed setting achieves. 
At high doppler frequency, the compared algorithm may be 
worse than the fixed setting. 

The performances of these three algorithms when the doppler 
frequency varies with time are shown in Table I Villi The doppler 
frequency changes every KM frames and it varies from 2 Hz, 
4Hz, ... to 20 Hz randomly. 

C. Analytical Results 

Here we provide the analytical results for comparison. Figure 
[Tol shows the analytical value of PRADA performance. It is 
observed that these curve match that in figure [9] Table [IX] 
also show the analytical value of PRADA performance. The 
throughputs when doppler frequency is 4Hz match that in Table 

VI. Conclusion 

In this paper, the problem of the feedback frequency reduction 
has been formally addressed and investigated. A close form 
expression of expectation FER over arbitrary M frames is 



■G — PRADA-A (M=120) 
■B— PRADA-B (M=30,K=4) 
V— QPSK (2/3) 




10 12 14 16 
Doppler Frequency (Hz) 



Fig. 10. Analytical results of throughput with doppler frequency variation. 



Doppler frequency (Hz) 


PRADA-A 


PRADA-B 


2 


1911.65 


1939.18 


4 


1831.16 


1872.55 


6 


1779.94 


1828.52 


8 


1745.16 


1798.56 


10 


1742.85 


1779.25 


12 


1737.48 


1764.22 


14 


1733.25 


1754.44 


16 


1730.02 


1747.02 


18 


1727.65 


1740.13 


20 


1725.66 


1736.69 



Throughput unit: data bits/frame 
TABLE DC 

Analytical results of throughput with doppler frequency 
variation. 



provided. On the foundation of the derivation, we were able 
to propose two adaptive algorithms, PRADA-A and PRADA-B 
with different degrees of feedback reduction. It was shown that 
the first algorithm proposed achieved significant performance 
improvement over fixed modulation. The second algorithm 
further reduced the frequency of feedback, moreover, introduced 
a novel adaptive transmission technique combing the feedback 
CSI and FER. The algorithms introduced in this work allow 
the system designer to flexibly switch the adaptation period 
according to the application under consideration. 



Appendix 



A. Derivation of {F™ k (s r , I)} 

First we define a random variable 



X, 



0, successful transmission 

1, erroneous transmission 



Recall that Ef* is the FER when using the setting s r in the 
channel state lUj. Given the used setting s r and the channel state 



(9) 



Wi, the probability generating function (PGF) of Xj\s r ,Wi is 
<j>x } \ Sr , Wi {u) = ^ x P[Xj=x] 

Denote the channel states in M frame transmissions as a random 
vector w = [w^w^ . . . wt M ]. The conditional PGF of the total 
number of erroneous frame in M frame transmissions is 

4>X 1 +X 2 + ---+X M \s r , V r{u) 

= 4>X j |8 P ,u> 4l \s r ,w i2 (w) • • • 4>X 3 \s r ,w iM M (10) 

We need to take expectation of the conditional PGF with respect 
to the probability of every possible w. This can be done by 
combining transition matrix and {ip SrtWi (uj)} as 



o 



where 



o 



















1 N,N 











(ii) 



(12) 



We can find the PGF of erroneous frame number in a setting 
adaptation period by matrix multiplication. That is, 



h na(u) H s N r 2 (uj) 



i,iv( w ) 



h£ t (u>) 



H 



N,N 



(13) 



Note that H?r(u) is simply the PGF of F^(s r , I). That is, 

M 

^M = EO S -0^- d4) 



1=0 



One can find {F^(s r ,l)} easily by H%(u) 
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